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ABSTRACT 



1. OBSERVATIONS 



X 



Comet Hale-Bopp was observed five times with 
ISOPHOT, the photometer on board the Infrared 
Space Observatory (ISO). Each time broadband pho- 
tometry was performed using 4 different detectors, 5 
apertures and, 10 filters covering the range between 
3.6 to 175 /Jin. Calibration measurements using the 
internal Fine Calibration Source were done together 
with each individual measurement. Background ob- 
servations were performed with identical instrument 
settings at the same positions on the sky several days 
after the comet observations. The observation strat- 
egy and the initial data reduction steps are described 
in some detail and the resulting in-band power values 
of the Hale-Bopp observations and their uncertainties 
are derived. Initial reduction of these measurements 
was performed in 3 steps: (1) processing of raw data 
by removing instrumental and energetic particle ef- 
fects, (2) averaging of the individual signals, and (3) 
determination of the detector responsivities and their 
uncertainties. The detector signal is determined by 
two different methods and the best value is chosen. 
At the present level of processing uncertainties range 
from 10% to a factor of 3 for the low power levels at 
short wavelengths. The in-band power levels at dif- 
ferent wavelengths varied over 3 orders of magnitude. 



*ISO is an ESA project with instruments funded by ESA 
Member States (especially the PI countries: France, Germany, 
the Netherlands and the United Kingdom) and with the par- 
ticipation of ISAS and NASA. 



Comet Hale-Bopp was a unique target-of-oppor- 
tunity comet during the lifetime of the ISO satellite. 
Photometric observations of comet Hale-Bopp were 
performed with the ISOPHOT instrument (Lemke ct 
al., 1996). The aim of this paper is to describe the 
observations and the initial data reduction scheme 
in enough detail in order to understand the present 
uncertainties of the measurement results. Astronom- 
ical interpretation of the results in terms of the spec- 
tral energy distribution requires additional color and 
aperture corrections and will be reported in a sepa- 
rate paper (Grim et al., in preparation). ISO viewing 
constraints together with the comet's orbital geome- 
try limited the visibility of Hale-Bopp for ISO. Ob- 
servations were only allowed in the solar elongation 
range of 60° to 120°. In addition, the Earth had to 
be at least 24° away from the satellite axis. Further 
pointing constraints were imposed by the Moon and 
Jupiter. Because of these constraints, comet Hale- 
Bopp was only observable for ISO during two peri- 
ods in spring and autumn 1996 and again in win- 
ter 1997. Hale-Bopp's perihelion passage occurred in 
April 1997. 

Within the visibility period of Hale-Bopp the low- 
est possible background had to be found in front of 
which the observations were to be performed. The 
search for the lowest background along the arc of 
Hale-Bopp's orbit was done in two subsequent steps. 
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First, a raw selection based on the IRAS all-sky sur- 
vey (ISSA) plates was performed. 'Dark regions' at 
12, 25, 60 and 100/im along the cometary track on 
the ISSA plates were selected. Second, the selected 
regions were examined with IRSKY (IPAC) to assure 
that no known point sources were within 5 arcmin of 
the cometary track, and no prominent (cirrus) struc- 
tures were visible on the 100/im ISSA plates. Also 
the background estimates for each filter of the ISO 
observations were derived with IRSKY. 

Background subtraction is crucial for photometric 
measurements but standard simultaneous offset mea- 
surements were dismissed because of the extended 
coma of Hale-Bopp. Instead we introduced so-called 
'shadow measurements'. A shadow measurement is 
the repetition of the Hale-Bopp measurement, called 
master observation, at the same position on the sky 
that was tracked in the master observation, at a time 
when the comet and more importantly the coma had 
moved away from that sky position. Shadow observa- 
tions were typically performed about one week after 
the master observations. For this kind of background 
measurements we assumed that the seasonal varia- 
tions of the zodiacal emission are small. Because 
of scheduling problems the 3rd Hale-Bopp observa- 
tion was not executed as planned but was delayed by 
about one month and, therefore, these master and 
shadow observations are a few degrees apart. 

Four different detectors of ISOPHOT (Lemke et al, 
1996) were used for the observations (Table |l|): de- 
tector PI, a Si:Ga-detector for observations up to 
16 /im, detector P2, a Si:B-detector for the 25 /im ob- 
servations, the C100 camera, a 3x3-pixel-array of nine 
Ge:Ga-detectors for 60 and 100 /im observations, and 
the C200 camera, a 2x2 array of four stressed Ge:Ga- 
detectors for the 175 /tm observations. All observa- 
tions were done in single pointing, single filter abso- 
lute photometry mode (AOT P03 and P22), i.e. with 
a subsequent measurement of the Fine Calibration 
Source (FCS). Multi-aperture measurements (AOT 
P04) of comet Halo-Bopp are reported in a compan- 
ion paper (Peschke et al., 1998). 



2. DATA PROCESSING 



Processing ISOPHOT data from raw data to mean- 
ingful astronomical values is a multi-step process 
with several decisions to be made that affect the end 
result. This process is described in some detail in 
order to get an estimate of the uncertainties of the 
result, and, to allow other researchers to compare 
their results with ours. For the initial data reduc- 
tion - that is described in this paper - the ISOPHOT 
Interactive Analysis tool (PIA, Gabriel ct al. 1997) 
is used as a baseline. In principal, it takes care of 
detector and instrumental effects and finally results 
in in-band power values for each individual comet 
and background measurement (per detector pixel and 
per bandpass). The second step that comprises color 
corrections, corrections for flux deficits due to inac- 
curate pointing and correct aperture scaling for the 
distributed coma brightness will be discussed later 
(Grim et al., in preparation). All these corrections 
are needed to derive a spectral energy distribution 
in a standardized aperture over a wide wavelength 
range. 



The five ISOPHOT Hale-Bopp observation sequences 
reported here comprise 558 individual measurements, 
counting all object, shadow, and calibration (FCS) 
measurements at up to 10 wavelengths, some of which 
were multi-pixel measurements with up to 9 pixels. 
Each of these measurements had to undergo all initial 
steps of data processing described below before they 
were combined into 44 Hale-Bopp and background 
signal values (Table [l]) . 

During the exposure to an infrared source the 
ISOPHOT detector output voltage ramps-up with 
time at a rate that is a function of the infrared flux 
received by the detector. This voltage is recorded 
and read-out in preset time steps of 1/128 s to 1/4 s. 
The output voltage is digitized with 12 bits resolution 
over a dynamic range of about 2.2 V. Normally, be- 
fore the detector output voltage reaches saturation, 
after a fixed number of read-outs, the voltage is reset 
and a new voltage ramp is started. Each individual 
measurement consisted of 4 to 683 voltage ramps. 
That way, from 512 to 4096 individual voltage val- 
ues per measurement were obtained. The individual 
read-out programs were defined in the observation 
planning phase dependent on the expected flux lev- 
els. Table [j] shows parameters of the different mea- 
surements. 

The objective of the initial data processing is to de- 
rive the flux dependent slope of the voltage ramps 
[V/s] and to determine the responsivity of the de- 
tector. The whole data was reduced homogeneously 
with PIA version 7.1. At raw data (ERD) level, non- 
linearities of the detector response and of the elec- 
tronic amplification are corrected by using the al- 
gorithms provided by PIA. During the whole data 
reduction process, the default criteria for discarding 
data points have been applied. In case of only 4 read- 
outs per ramp, just the two middle data points are 
left for the slope determination of the ramp. 

Next, the effects of high energetic particle hits are 
removed which cause the signal to display discrete 
steps of variable amplitudes on top of the linear rise. 
These glitches are removed in two steps, partly on 
the raw data level, partly after the determination 
of the slopes. The slopes of the integration ramps 
have been determined 1) by first order fits to the 
ramps (method 1) and , 2) by determining the in 
pairs-differences of subsequent non-destructive read- 
outs (method 2). By applying the drift recognition 
algorithm with its default settings, an average signal 
for the measurement was deduced. 

To get confidence in the values derived with these 
two different methods within PIA, a third indepen- 
dent method was introduced to have an independent 
handle on the effect of glitch removal as well as on 
detector drifts. For this, only the non-linearity cor- 
rection was applied to the raw data and further pro- 
cessing occurred outside PIA. All in pairs-differences 
were computed, regardless of the type of readout. At 
this stage, the overall trend in the data, the effect 
of detector drifts and glitches was clearly visible in 
the data. Glitches could easily be identified by their 
distance to the bulk data, based on very good statis- 
tics. Next, a myriad algorithm was applied to the 
data which computed and analyzed the distribution 
of slope values. A trim fraction was introduced to re- 
move the tails of the distribution that contain mainly 
glitches. Then the remaining data is analyzed for its 
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momentum and a stable mean method is applied. 
The derived value reflects the average signal of the 
measurement. This value together with the appropri- 
ate error is written back into a structure for further 
processing within PIA. 

Method 3 is superior to the previous methods when 
the measurement comprises only a small number of 
ramps. We calculate the slopes by all methods in 
parallel but use method 3 for the final result when- 
ever possible. However, in case the voltage increase 
per read-out is less than 3 mV method 3 deteriorates 
because of digitization noise and results from method 
1 or 2 are used. Independent of the method by which 
the average signal slope was derived, the standard 
reset interval correction is applied, dark current is 
subtracted and the effect of vignetting is corrected in 
case of the C detectors. 

All the above processing steps are applied both to 
the object and the subsequent FCS measurements. 
In an FCS measurement the detector is irradiated by 
the internal calibration lamp (FCS) and the actual 
detector responsivity can be determined. However, 
this method works only reliably when the infrared 
power onto the detector from the FCS is within about 
a factor of 3 of the power from the object. The 
FCS power was set beforehand on the basis of es- 
timates of the object brightness. In case of larger 
deviations - this was the case in about 2/3 of the 
measurements - the default responsivity should give 
more reliable results. The actual and default respon- 
sivities sometimes showed large deviations: up to a 
factor of 3 and more, especially in the case of the very 
low power values at short wavelengths. In several of 
these cases the FCS heating power and the corre- 
sponding power of infrared radiation were out of the 
calibrated range. In any case, we use the deviation 
of the different responsivity values (default, and ac- 
tual responsivities determined by all three methods) 
as a measure of the uncertainty of the responsivity 
and the derived in-band power value. Figure nj shows 
the derived in-band power values of the Hale-Bopp 
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measurements. For a single observation the in-band 
powers at different wavelengths vary over about 3 
orders of magnitudes. Note especially the close simi- 
larity of the in-band power values of the 27-SEP-1996 
and 7-OCT-1996 observations, except for the value at 
175/im. Due to scheduling problems the positions of 
the master and the shadow observations were offset 
by a few degrees which only showed a strong effect 
at the longest wavelength. 
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Figure 1. In-band power values of the five Hale-Bopp 
observations. Different symbols characterise different 
observations: +: 25-MAR-1996, *: 27-APR-1996, o: 
27-SEP-1996, A; 7-OCT-1996, □; 30-DEC-1997. 
Two measurements at 25fim wavelengths were not 
completed because of spacecraft problems. 
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Table 1. Log of Hale-Bopp and corresponding shadow observations. For each individual ISOPHOT measurement, 
detector, wavelength of filter, aperture size and integration time are given. For some observations the complete 
set of filters could not be used. For the specific measurements (Hale-Bopp and shadow) the number of voltage 
ramps, number of read- outs per ramp and the resulting responsivity and in-band power values are given 
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